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Abstract 

Superconductivity in Ceo materials is modeled in terms of an in- 
tramolecular, nonretarded attraction. It was shown previously that, at 
intermediate coupling, the model possesees a state in which s-wave super- 
conductivity coexists with a charge density wave, the latter stabilized by 
intersite repulsion. The CDW causes Tc to decrease near half filling. We 
argue that hole-doped Ceo, for which Tc peaks away from half filling, is a 
possible candidate for this state. But the electron-doped Ceo and AaCeo 
are conventional superconductors, stabilized against a CDW by metallic 
screening, which is treated in a parameter-free fashion. 

PACS: 74.10 +v, 74.20 -z, 74.70 Wz, 74.72 -h 

About a decade ago, superconductivity with surprisingly high Tc's (30— 40A') 
was observed in A^Ceo, where A is an alkali atom p 
TcS {52 K — 117K) have been seen in hole-doped Cgo 

Cooper-Schrieffer (BCS) picture, origin of such high Tc's has been attributed to 
special properties Cgo, such as, high density of states (DOS) and large phonon 
frequency. However, variation of Tc with the density n does not follow the DOS. 
Also, in the electron-doped Cgo and AaCgo the maximum of Tc occurs near half 
filling. But in the hole-doped material, Tc peaks at a density of rt = 3 — 3.5 
holes/molecule, and decreases rapidly toward half filling (5 holes/molecule) 
Such a behavior has been predicted previously on theoretical grounds H, 1). The 
sharp decline in Tc is caused by the appearance of a charge density wave (CDW) 
near half filling which, due to the frustration effects of the lattice, can coexist 
with superconductivity. Here we examine the conditions for the existence of 
such a state in Cgo- 

The possibility of a CDW instability in Cgo was recognized early by Zhang, 
Ogata and Rice|^. As discussed below, it becomes more likely as Tc increses. 
The key issue is the effect of Coulomb repulsion which opposes superconduc- 
tivity, but favors the CDW. We will focus mostly on the normal state, and 
consider realistic potentials. Wc find that if there is no metallic screening, the 
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energy gained by the CDW increases monotonically with increasing interaction 
strength. However, in the presence of metalHc screening, the gain has a maxi- 
mum. It is too small to stabilize the CDW state in the lower Tc materials such 
as AaCgo- But for the hole-doped case, the CDW state can not be ruled out. 
We discuss some normal-state anomalies of this state. 

Crystalline Cgo has a face-centered-cubic (FCC) structure, and is a semicon- 
ductor. The conduction band is formed from the 3-fold degenerate ti^ molec- 
ular orbitals, and the valence band from the 5-fold degenerate orbitals. In 
A3 Ceo, each alkali atom donates one electron to the conduction band making 
it half-filled. Large molecular size leads to a large lattice spacing (a ^ 10 A), 
and consequently very narrow bands of width W ^ 0.5eV, which can thus 
be adequately described by a tight-binding model with nearest-neighbor hop- 
ping. This implies that geometric effects such as frustration and nesting are 
important. 

Our theory is based on two key assumptions. (1) Low energy physics of 
the doped system can be described by an effective Hamiltonian defined entirely 
within the relevant band, i.e., the hiu band for hole-doped Cgo, and the tiu 
band for electron-doped Cgo and A^Cqq. The effect of integrating out other 
bands and phonon degrees of freedom is to renormalize the interactions. (2) 
The attractive interaction responsible for superconductivity is intramolecular in 
origin and effectively nonretarded, i.e, the characteristic frequency loq is compa- 
rable to the Fermi energy. Indeed, intramolecular phonons that are thought to 
be responsible |^ have rather high frequencies: htoo ~ O.lbeV — 0.2eV, easily 
comparable with the Fermi energy. 

The model always has a solution in which the normal state is conventional, 
and superconductivity is BCS like. In the BCS theory Tc is given by 

kTc w l.Snwoe"^/^, (1) 

where A ~ U piep) is the dimensionless coupling constant, p{tF) is the density 
of states (DOS) at the Fermi level, U is the attractive interaction. The observed 
high Tc's are usuallyl^ attributed to: (1) a relatively high DOS due to the nar- 
row bandwidth (p ex W~^) and orbital degeneracy, and (2) a large prefactor due 
to high phonon frequency. Furthermore, in AaCgo, Tc increases with increasing 
size of the alkali atoms. This is because as the lattice expands, W decreases. 
Similarly, the five-fold degeneracy of the hiu level should lead to higher DOS 
, which will account for the higher in the hole-doped material. 
However, not everything fits. According to the simple BCS formula, the 
n-dependence of should track that of the DOS (through the n-dependence of 
the Fermi energy). This is actually not the case. For the FCC lattice the DOS 
has sharp features. The measured Tc is smoother, its maximum is far removed 
from that of the DOS, both for the electron doped or hole doped materials]^. 

We point out that the problem is partly due to the assumption that the 
attractive interaction is retarded, i.e., confined within Tiloq << ep of the Fermi 
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surface. Then one can approximate the DOS by piep), and use hujo as a cut-off 
in the energy integral which accounts for the prefactor. This procedure is not 
vaUd when fkoo is comparable to e_F. In fact, as shown previously in the 
nonretarded case Tc varies much more smoothly with n. The reason is that for 
general values of huJo , a more accurate approximation is given by 

where hwi = min{eF, fn-Oo), hw2 = min{W — ep, fn-Oo), and ep is measured from 
the bottom of the band. Here po is the weighted average of the DOS p{e) over the 
allowed region with a slowly varying weight factor tanh(e/2fcrc) (normalized 
to unity). For small fnuo we recover the usual BCS formula, whereas for large 
huJo, the variation of Tc with n is smoothed out considerably. In this 
better estimate is po ~ 1/W. 

A more important issue is the possibility of a CDW instability which can 
not be ignored when Kloq is comparable to ep. We consider a one-band model, 
characterized by a single-particle energy et, an on-site attraction [/, and an 
intersite repulsive interaction V{r). The relevant dimensionless energy scales 
are Up and Vp. The effect of orbital degeneracy can be approximately taken 
into account by changing p, or, equivalently, by changing U /W and V/W. The 
possibility of a CDW was studied in reference using this model, but without 
the V term. However, as shown previously, inclusion of even a nearest neighbor 
V leads to nontrivial consequences ^ . We now generalize to V{r) of arbitrary 
range. 

In the simplest mean-field approximation, and in the absence of supercon- 
ductivity, a CDW with a wavevector Q would appear for 

{U + 2Vq)x{Q)>1, (2) 

where x(Q) is the charge susceptibility for the noninteracting system, and Vq = 
— ^j._^Q e*'^ '"F(r). The minus sign is included in order to make Vq positive. 
For a bipartite (e.g., simple cubic, BCC etc) lattice, there is a Q for which 

Ek+Q = -ek (3) 

for all k. Then the Fermi surface is perfectly nested at half filling, leading to a 
logarithmically divergent x at T = and a CDW instability for any {U + 2Vq) > 
0. For the nonbipartite FCC lattice, x(Q) is finite even at half filling as there 
is no nesting, i.e., no Q for which Eq.(3) holds. A conventional metallic state 
is then stable for {U + 2Vq) < l/x(Q) ~ 0.375W. Above this, a CDW state 
appears with Q = (0,0, 27r/a), which corresponds to a CDW with a density 
n(r) = n + nie'*^ "", i.e., on successive planes perpendicular to the z axis, the 
density alternates between n + ni and n — ni. 

From (2), it appears that a CDW can exist even V — 0. However, such 
a state is unstable against s-wave superconductivity, since the pair suscepti- 
bility is logarithmically divergent. The inclusion of V{r) changes this picture. 
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due to two reasons. First, the V term favors the CDW energetically. At the 
mean-field level, the excess potential energy per site in the CDW state equals 
— {U + 2VQ)nl/4:. Physically, the CDW forms an "ionic crystal" with Vq propor- 
tional to the appropriately generalized Madelung energy However, compared 
with a bipartite (e.g., BCC) lattice, charge ordering on an FCC lattice is frus- 
trated, leading to a somewhat lower energy gain. Nonetheless, we have verified 
numerically that Vg > 0, for the screened Coulomb potentials (Fig 1) which are 
the potentials of interest. Second, V{r) always opposes superconductivity by 
keeping particles apart, stabilizing the CDW state at larger Vq. 

The CDW splits the noninteracting band into two bands of energy 

£;±(k) = ei(k)±[e^(k)+A^]^ (4) 

where ei(k) — ^(ek + Ck+q) and e2(k) — ^(ek — £k+q)- At the same k, the 
bands are separated by a gap Ac = ni{U + 2Vq)/2. For a bipartite lattice 
the geometrical property, Eq. (3), ensures that ei(k) — 0, leading to insulating 
behavior. However, Eq.(3) does not hold for the FCC lattice. Then for small 
enough Ac, the maximum of the lower band lies above the minimum of the 
upper band which, in general, is at a different k point. This leads to a semimetal 
with a partially gapped Fermi surface. The is a direct consequence of geometric 
properties frustration and (lack of) nesting. Eventually, there is a transition to 
a an insulating state for {U + 2Vq) > 0.562W^. 

Furthermore, residual attractive interaction within the semimetal leads to 
superconductivity. For a nearest-neighbor V, the CDW has been shown to 
coexist with superconductivity over a substantial region in the parameter space 
The presence of the CDW lowers T^, and yields |^ a T'c(n) which is very 
different from that in the conventional state (no CDW). In either case, Tc does 
not track the DOS. In the conventional state, it peaks near half filling and 
gradually decreases to zero at the band edges. This is similar to the behavior 
observed in electron doped Ceo and A3 Ceo- 

By contrast, in the CDW-superconductor Tc has a minimum near half filling 
where the CDW instability is strongest. Away from half filling, Tc increases 
as the CDW weakens. Eventually, there is a transition to the conventional 
state. Consequently Tc reaches a maximum near the transition point and then 
decreases following the behavior in the conventional state. This is very similar 
to the behavior observed in the hole doped Cgo- So far only one maximum has 
been seen because experiments have been limited to only one side of half filling. 
The theory predicts two, one on each side. These results should be qualitatively 
correct for the more realistic screened Coulomb potential, although one expects 
a more rapid decline of Tc toward the band edges since screening length increases 
with distance between particles. 

Are the observed Tc's consistent with the theory? Consider the conven- 
tional state at half filhng which exists for {U + 2Vq)/W < 0.375. To esti- 
mate U /W from the measured Tc, we use formula (1), conveniently rewritten as 



4 



1/A = —ln{kTc/1.3huj(j), to compute A, which is treated as a phenomenological 
coupling constant. At half filling, we can take t^o to be the (intramolecular) 
phonon frequency, which is same for all materials. Noting that A varies slowly 
with luq, we set huo = ISOO-ftT, which is reasonable. For the electron doped Cgo 
(Tc UK) this gives A = 0.146. However, for RbaCeo {Tc = 33K), A = 0.245. 
Next step is to relate A with the theoretical parameter U /W. As discussed 
before, in the absence of ^(r), and for large hujQ, A ~ U/W. However, when 
V'(r) is included, A will be reduced so that U/W > A. This means that electron- 
doped Cgo and A3C60 can be conventional metals, but for the larger alkali (Rb 
etc) materials Vq/W has to be rather small. 

For the hole-doped Cgo, the maximum Tc of 52K yields a larger value of A 
(« 0.28) which could be sufficient for the CDW state. However, the maximum 
is away from half filling which in our theory is due to the transition to the 
CDW-superconducting state. In order to be consistent, we have to use the Tc 
of the corresponding conventional state at half filling. This is about 30 to 50 
percent higher, which gives U/W > X ^ 0.32 — 0.34. For any reasonable values 
of Vq these materials should be CDW metals near half filling. Any expansion 
of lattice will increase U/W, making the CDW even more likely Q]. 

The intriguing question is, why is A3C60 not a CDW metal, as suggested 
previously^? For example, for RbaCgo to be conventional, Vq/W has to be 
~ 0.05 or less, almost two orders of magnitude smaller than the bare Coulomb 
energy scale e^/a « 2.6W. Now, there are two sources of screening. (1) Screen- 
ing by other molecular bands which can be approximated by a static dielectric 
constant e. One expects e to be somewhat larger than that in ordinary semi- 
conductors because of small (molecular) band gaps. (2) Since the transition is 
from one metallic state to another, we must include metallic screening which 
leads to exponentially short-range potential and much smaller Vq. To see this, 
we consider the screened Coulomb potential 

V{r) = -e--/^ (5) 
er 

where ^ is the Thomas-Fermi screening length. 

Let us define the dimensionless quantity a by writing Vq — voct{(./a), where 
Vq = e^/ea sets the energy scale for the long-range part. Then a{x) is generalized 
Madelung constant [||, defined for screened-Coulomb potential. As shown in 
Fig. (1) a, which measures the energy gained by the CDW, increases with 
the range of interaction ^/a. In the absence of metallic screening, ^ = 00, 
and a ~ 1.6. With decreasing ^, a decreases rapidly below ^ ~ a, which is 
essentially the metallic region. Fig. 1 also shows that a is always smaller for 
the FCC lattice than the bipartite BCC lattice. This is due to frustration which 
is another — albeit smaller - source of reduction of Vq . 

It is possible to have a CDW state in every case by choosing the param- 
eters Vq (i.e., e) and ^ appropriately. However, a much better parameter-free 
understanding of screening can be obtained as follows. First, note that ^ itself 
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depends on vq . The Thomas-Fermi expression for ^ [ |lO[ , written in the lattice 
language, reads: 

^/a=[{in^vo)^r^/^ (6) 

where fj, is the chemical potential, and the geometrical parameter 7 — \/2 for 
the FCC lattice. At low T, dn/dji w pi^p) ~ 2/W. Let us define an effective 
bandwidth by 2/Weff — dn/dp,. Since a depends only on ^/a, we see that 
the dimensionless quantity Vq/M^e// is a universal function of the dimension- 
less coupling constant vg/Weff only. This is plotted in Fig. 2. Notice that 
Vq/We// ~ Vq/W, the quantity of interest. Naively one would expect Vq to 
keep increasing as vq increases. On the contrary, Vq/VFe// has a broad maxi- 
mum. This is because an increased vq (i.e., a larger charge) makes the screening 
length ^ shorter (Eq. 6) which decreases a. The maximum occurs at pvo ~ 0.15 
and has a value of about 0.042. This is a pure number, which depends only on 
the lattice type, but not any parameters. Hence, we have the remarkable result 
that irrespective of how large or small vq (i.e., e) is, 2Vci/Weff < 0.084. The 
smallness of this number, combined with the above estimates of U/W, explains 
why A3 Ceo is conventional. We stress that this drastic reduction of Vq would 
not occur if the CDW state were an insulator since there is no metallic screening 
at the transition. 

Since Vq/VFe// varies slowly (fig. 2), its maximum value also provides an 
estimate for its magnitude. For example, 2VQ/Wi,ff > 0.06 in the entire range 
of 0.05 < Vo/Weff < 0.4 (52 > e > 6.5). This is sufhcient for the hole-doped Cgo 
to be a CDW semimetal. This analysis is of course approximate. One way to 
see if the state is experimentally realized is to look for normal state anomalies. 

The semimetal is described in terms of two bands separated by a direct gap 
Ac{T), as well as a Fermi surface. The gap vanishes at Tcdw, above which 
the normal state is conventional. Below Tcdw, the band parameters become 
temperature dependent through Ac{T), leading to distinct features. Thus, the 
carrier density becreases with increasing Ac, leading to peak or a shoulder in 
the resistivity as a function of T. Other quantities, such as the specific heat, also 
show similar behavior. To see this, consider the approximate low T expression 
for energy 

E{T)^Eo{Ac) + ^{kTrp{eF), (7) 

where Eq(Ac) is the ground-state energy and is the corresponding Fermi en- 
ergy for a given Ac{T). For T > TcdWi Ac = 0, and Eq and p are independent 
of T, leading to linear heat capacity and T-independent magnetic susceptibility. 
Below Tcdw, Eq and p(ep') become temperature dependent through Ac(T). 
We have solved the mean-field equations numerically at finite T. Fig. 3 shows 
typical behavior of specific heat as a function of T /Tcdw, which is linear in T 
above T > Tcdw, and again becomes hnear at low T, as Ac(T) saturates, but 
with a smaller slope, corresponding to reduced value of p{£f)- The sharp feature 
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below TcDW is due to the rapid decline in p{eF) as function of temperature, 
also shown in Fig. 3. These features should be experimentally observable in the 
region where Tcdw > Tc, i.e., close to half filling. 

In conclusion, we have shown that despite relatively high Tg's, frustration 
and metallic screening suppresses a CDW instability in A3 Ceo- But for hole- 
doped Ceo, a mixed state can not be ruled out. If the anomalies discussed above 
are not observed, then the CDW is absent. This is possible, given the simplicity 
of the model and approximate nature of the analysis. Even then the instability 
can not be be very far away, and can occur as Tc is increased further by lattice 
expansion. The author thanks C. Jayaprakash and T. L. Ho for discussions. 
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FIGURE CAPTIONS 
Fig.l: The generalized Madelung constant a = Vq/vo = YIt^o '' (r/a) — plotted 
against dimensionless screening length ^/a, where a is the nearest neighbor 
distance. Note that a depends only on ^/a and the lattice type. The solid line 
is for the FCC lattice and the dashed line for the bipartite BCC lattice, showing 
that energy gain from the CDW is less for the FCC lattice due to frustration. 
In usual metals, ^/o < 1 . 

Fig. 2: The dimensionless quantity VQ/W^ff is plotted (solid line) against 
vo/Weff where We// = 2/^ ~ W. Also shown (dashed line) is 10^/a which, 
like a, is slowly varying except for very small Vq. 
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Fig. 3: Scaled specific heat (solid line) and density of states at the Fermi level 
(dashed line) are plotted vs T/Tcdw for the CDW- metal at half filling. The 

parameter values are t/+2V^ = 6.71i, Tc = 0.768t, The prominent feature below 
TcDW is due to the CDW and reflects the rapid rise in the DOS. 
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